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Abstract 

I discuss in these lectures vortex-like classical solutions to the equa- 
tions of motion of gauge theories with spontaneous symmetry break- 
ing. Starting from the Nielsen-Olesen ansatz for the Abelian Higgs 
model, extensions to the case in which gauge dynamics is governed 
by Yang-Mills and Chern-Simons actions are presented. The case of 
semilocal vortices and also the coupling to axions is analyzed. Finally, 
the connection between supersymmetry and the existence of first order 
BPS equations in such models is described. 
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Turning and turning in the widening gyre 
The falcon cannot hear the falconer 
W.B.Yeats, The Second Coming[T] 

These notes cover the topics of 4 lectures on vortex solutions in 
spontaneously broken gauge theories and their connection with super- 
symmetry, given at the 4th Chilean School of Astrophysics, Cosmology 
and Gravitation held in Valparaiso last October. The bibliography at 
the end of the notes is not complete: it only refers to works containing 
results discussed in the lectures. A more complete list of references 
can be found in some of the excellent books and reviews on this topic 

1 Nielsen Olesen vortices 

Vortices are very ubiquitous objets. Starting from the 1957 work of 
Abrikosov [8j on superconductivity, belonging to the area of condensed 
matter, they were then rediscovered by high energy physicists in a 
field theory context, after the well-honnored paper by Nielsen and 
Olesen [9]. In the present century a revival took place due to their 
cosmological interest in connection with cosmic superstrings (see [10] 
and references therein). 

Nielsen and Olesen [9] introduced vortices in high-energy physics 
in an attempt to describe classical string-like structures with the prop- 
erty of having equations of motion identical with those of the Nambu 
dual string. 

In their work, Nielsen and Olesen started from a 3 + 1 Lagrangian 
for an Abelian U{1) gauge theory coupled to a charged scalar with 
a symmetry breaking potential and first found a classical string-like 
(vortex) solution. Then they obtained the remarkable result the action 
is proportional to the area of the surface swept out in space and time 
by such vortex. This result is the basis of their identification of the 
Abelian Higgs Lagrangian with the Nambu Lagrangian for the dual 
string [TT] . 

By making a static axially symmetric ansatz they reduced the 
Euler-Lagrange equations to a pair of second order non-linear coupled 
differential equations in the radial variable. Although they didn't find 
exact solutions (in fact, analytical vortex solutions are not known) 
they gave arguments to show that nontrivial solutions exist. They 
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correspond to vortex-like object carrying a magnetic field in its core, 
where the scalar field vanishes. The magnetic fiux is quantized (al- 
ready at the classical level). 

Such solutions look like type-II superconductivity vortices and this 
is no coincidence: it is not difficult to see the connection between the 
Abelian-Higgs model Lagrangian and the Ginsburg-Landau fenomeno- 
logical model for superconductivity, provided the Higgs scalar in the 
relativistic field theory is identified with the order parameter in the 
Ginzburg-Landau model. In order to have type-II superconductivity 
the relation between the gauge coupling constant and the symmetry 
breaking coupling constant A should be < 8A (with the conventions 
adopted below); for > 8A a complete Meisner effect takes place 
(type-I superconductivity). The limiting point 

= 8A (1) 

plays also an important role in the field theory /string theory context. 
In fact, a surprising connection was already notice in 1976, when it 
was observed pL2j that relation ([1]) is precisely the same needed for the 
supersymmetric extension of the model |13j-|14j. 

The Abelian Higgs model dynamics in (3-|-l)-Minkowski space is 
governed by the action: 

S = Jd'x (^-\f,,F^'^ + l{D,^nD^<p) - VM)^ (2) 

Here F^^i, is the field strength, 



(p is a complex scalar. 



vm) = \{\<t^\''-<g) 



(3) 

(4) 

(5) 
(6) 



the covariant derivative is defined as 



and the symmetry breaking potential taken in the form 



Here (pQ is a real positive constant with [mass] dimensions. 
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The resulting equations of motion are 




with the Higgs field current given by 

ju = '-^{c^d^r - cj)*d,cP) + e2|</.|2yl, (8) 

The Nielsen-Olesen strategy to construct nontrivial regular solu- 
tions to the equations of motion with finite energy (per unit length) 
starts from a trivial (constant) solution and implies the following steps: 



1. Trivial zero-energy solution \(p\ = (f>Q , Ai = 



2. non-trivial but singular solution (fluxon with N units of magnetic 
flux) which in polar coordinates reads 

4> = (t>o exp{iNip) , Ai = Ndup , SijFij = ^■kN6'''^\x) 

^jeijFij = NeZ (9) 

This configuration is a gauge transformed of the trivial one. The 
appropriate gauge group element to consider is 

gN{f) = exp{iNip) (10) 

which is ill-defined at the origin. 

3. Regular Nielsen-Olesen vortex ansatz 

(t> = f{r)gN{^) , Ai = a{r)^g]:l{ip)digN{^) (H) 

with 

/(0) = a(0)=0, /(oo)=0o, a(oo) = l (12) 

This ansatz leads to the same magnetic flux as that in ([9]) but there 
is no singularity in view of the boundary conditions that one imposes 
at the origin on / and a. Conditions at infinity guarantee that each 
term in the Lagrangian goes to zero at infinity. 
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Being axially symmetric, ansatz (jlip can be interpreted as a mag- 
netic tube with quantized flux (a vortex). Being z- independent, such 
a solution has infinite energy. However, if we define an energy E per 
unit length, conditions (fT2|) ensure its finiteness, 

E = Jd'x + 1 |^^^|2 ^ ^(|^|2 _ ^^2)2^ (13) 

Being z independent, such a static vortex solution can be also inter- 
preted as a soliton solution in 2 + 1 dimensions. Because it depends 
only on x, y, it can also be thought as an instanton solution in d=2 
Euclidean dimensions. In any two of this cases, space at infinity can be 
seen as a circle 5^. But also the gauge group, U{1), can be identified 
with a circle so that one can see 5 at as establishing a map 

gN-.Sl^ (14) 

Such maps are classified into homotopy classes Hi (where the subindex 
refers to the dimension n of the sphere at infinity, in this case n = 1). 
Now, one has 

Ui{S^)=Z (15) 

and in this way one connects the vortex magnetic flux Qwith a topo- 
logical charge. 

Unlike monopoles or instantons, no analytic solution to the equa- 
tions of motion is known. However, numerical solutions can be easily 
constructed and there are also proofs that solutions exist (see |3j and 
references therein). The asymptotic behavior of the vortex solutions 
is 

a(r) = \- mrKi(—mvr) 

e 

(j){r) = 00 + 0{exp{-m^r)) (16) 

where mv and are the masses of the vector and scalar particles 
respectively, 

mv = ecpo (17) 

m0 = 2V2X(t)o (18) 

Note that when condition ([T]) is satisfied both masses coincide, this 
being the first indication that the model can be viewed as describing 
the bosonic sector of an = 2 supersymmetric theory, with the scalar 
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and the gauge field in the same supermultiplet. In fact when my ~ 
we have a weh defined vortex fine or a weh defined string of extension 
l/my ~ Inside the string one basically has magnetic field while 

the scalar field vanishes. Outside the string it is the magnetic field 
which is nearly zero while the Higgs scalar is practically in its vacuum 
value 101 = 00- 

Explicit solutions in the whole r range can be more easily found 
observing that, under magic condition ([T]) the second order Euler- 
Lagrange equations can be reduced to first-order ones. There are 
different ways of arriving to this result. One was discovered in [12j 
and starts by demanding that the space-space energy-momentum com- 
ponents {Trr and T^^) vanish. The other one establishes a bound 
(Bogomol'nyi bound) for the vortex energy (per unit length) which is 
saturated whenever ([T]) holds [15]. Being configurations static, such 
a bound also corresponds to a bound for the action and hence they 
satisfy the equations of motion. 

The Bogomol'nyi bound 

Let us see how Bogomol'nyi found in ref.[T5] a bound for the energy 
starting from formula ([T|). To this end, let us write the scalar field in 
terms of its real and imaginary components 

c^={<p\4>'') = {<p-) , a = 1,2 (19) 
and scale fields and coordinates according to 

r^<l)Q(t)\ A,^<j)QAi, x^-^y (20) 

600 

SO that the energy reads now 




where the covariant derivative is defined as 

Df = J^^^i + e'^^Ai (22) 
One easily sees that the energy can be rewritten as 



6 



1 



2 

^\ (^^^JEiJ T ^e'''ei,di{rDj<p'') ] ] (23) 
1 



The term 



^e'^'e.AirDj^') (24) 

gives no contribution to the energy, as can be seen by using the Gauss- 
Stokes theorem and the fact that Dicf)'' = on S^^. Then, using ([9]) 
we have 

+ \ {D.,r T e'^'e^,!),/)' + - - 1)') ± '/'ovrA^ 

(25) 

Were the third term absent, the energy would just be the sum of 
to positive definite terms plus a term proportional to the topological 
charge. This is precisely what happens for the Bogomol'nyi point ([1]) 
for which one then has 

E > (PIttIN] (26) 

The bound is attained when the two perfect squares vanish, and this 
leads to the so-called Bogomol'nyi or BPS equations for vortices, 



Fij = TSij (0^ - 1) 



= ±e''^eijDj<l)^ (27) 



The "PS" refers to M. K. Prasad and C. M. Sommerfield who 
found exact monopole solutions for the non-Abelian version of the 
model for a particular (Bogomol'nyi) value of the coupling constant 
[16j (A = 0). That their solution corresponds to a Bogomol'nyi bound 
was first noted in [17] . 

As mentioned above, eqs. ()27p were independently obtained by ask- 
ing the stress tensor components Tij to vanish. This is consistent 
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with the eqs. of motion only if A = 8e^. When matter interaction 
exceeds that of the electromagnetic one (A > 8e^) vortices attract 
whereas in the opposite case they repel. In the former case a system 
of many superimposed vortices decays into separated vortices. In the 
case A = 8e^ vortices do not interact and this is consistent with thew 
fact that the energy is proportional to the topological charge. An ex- 
act (numerical) solution to the system (j27p was presented in [l2,j . The 
existence theorem for vortex (and monopole) solutions can be found 
in A. Jaffe and C. H. Taubes book on Vortices and Monopoles [3]. 

Chern-Simons vortices 

As explained above, if one defines the Maxwell-Higgs system in d = 
2 + 1 dimensions, vortices can be seen as static finite energy (solitonic) 
solutions of the equations of motion. Now, as it is well-known, there 
is in three dimensional space time a second Lorentz invariant term 
which can be added to the Maxwell term, namely the Chern-Simons 
term which, for a U{1) gauge theory reads 

Lcs = ^e'^^^F^^Ap (28) 

Here is a parameter with mass dimensions usually called the topo- 
logical mass. This name is due to the following important fact related 
to the Chern-Simons term [18j-|2U|. Consider the U(l) gauge theory 
with Lagrangian 

Lm+cs = -If^.F^^'' + l^e-'^^F^^Af, (29) 
Now, from the eqs. of motion for this Lagrangian, 

d^F'"' + fi*F'' = (30) 
with the dual field strength *F'^ given by 

= ^e-^PF^f, (31) 

one can find 

showing that the gauge field is massive. 
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Clearly the field equations are gauge invariant but the Lagrangian 
is not. Under a change 

A^^A^ + d^A (33) 

the Chern-Simons Lagrangian changes by a total derivative 

Lcs^Lcs + ^daCF^A) (34) 

The action is however gauge-invariant since the field strength fall off 
at large distances and one accepts appropriate A's. 

Let us then consider the Maxwell-Higgs model with the addition 
of a CS term, with Lagrangian 

S = Jd'x (^-\f^.F>^'^ + ^s^'^fF^^Af, + i(I)^0)*(D'^0) - VM) 

(35) 

The equations of motion now read 

flfF^^ + ^'F^ = i„ (36) 
D-A,* = ^ (37) 

Now, integrating the temporal component of eqs. ([5S]) (the Gauss law) 
in the plane one has 

J cfx {O'Fm + ne^JF'^) = J d^xjo (38) 

The first term in the l.h.s. is a surface term and gives no contribution. 
The second one is just the magnetic flux and that in the r.h.s. is 
identified with the electric charge of the configuration. One then has 

fi^ = Q (39) 

with 

If one looks for static configurations, the current component jo is just 

jo = e^\(l)fAo (41) 



9 



We then see that either 7^ or we cannot have vortex-hke solutions 
hke those discussed for the pure Maxwell-Higgs system. Interestingly 
enough, if vortex solutions exist, magnetic flux and electric charge 
are related according to (j39p . Being the magnetic flux related to the 
topological charge, the electric charge is then also quantized already 
at the classical level. 

In fact, it is very easy to extend the ansatz (jlip in order to include 
an ^0 [21], 

<f) = f{r) exp(m(/?) , Ai = na{r)dnp , ^0 = ao{r) (42) 

Apart from conditions (jl2p . one should ask 

00(0) = ao(oo) = (43) 

As it is to be expected in view of our previous discussion, the vector 
mass is modified due to the presence of the Chern-Simons term (while 
the scalar field mass remains unchanged) 

mA± = y?7^±| (44) 

= 2\/2A(/)o (45) 

Using the axially symmetric ansatz, numerical solutions can be found. 
The behavior of the Higgs scalar and the spacial components of the 
gauge field is qualitatively the same as in the pure Maxwell-Higgs 
case. Concerning the electric field, it vanishes at the origin and has 
its maximum value at finite r. The energy E of the vortex can be seen 
to be E± ~ log (e^/m^_|_) and hence the rriA- solution has lower 
energy. 

No Bogomol'nyi bound for the energy can be constructed when the 
CS term is included. Hence, there are in this case (and with this field 
content) no BPS first order equations for the Maxwell-Higgs-Chern- 
Simons model. 

There is however a way out to this problem; one is to modify the 
field content. The other one has been discovered in two simultaneous 
works j22j-|23j. Basically, it consist in not including the Maxwell 
term -that is, the gauge field dynamics is solely governed just the 
Chern-Simons action- and changing the symmetry breaking potential 
from the usual fourth order one to a sixth order potential which is 
acceptable from the point of view of renormalizability since the model 
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is defined in 2 + 1 space-time. The model can be seen as the truncation 
at large distances and low energies where the lower derivative Chern- 
Simons term dominates the higher derivative Maxwell term. 
The Lagrangian of the model is 



S 



Although one can include in the potential an arbitrary sum of fourth- 
order and sixth order potential, Bogomol'nyi equations will arise for 
a particular combination that can be written in the form 

^(') = ^l'AP(l0p-0g)' (47) 

There is no independent coupling constant in this potential because it 
has been already defined at the Bogomol'nyi point. Note that being 
defined in a d = 3 dimensional theory, it corresponds to a renormaliz- 
able potential. 

Variation of the action yield to the field equations 



e^f^Fp^ = f (48) 



2 

D^D^^ = -—— (49) 



where j" is still given by The time component of eq. (|48p is the 
Gauss law, 

/^i^i2 = JO (50) 

Noting that there is no metric in the Chern-Simons action (in this 
sense it is a topological Scwartz like action |24j ) one concludes that it 
does not contribute to the energy momentum tensor, which is defined 
as 

T.. = 2^ (51) 

Thus, the energy functional for static field configurations reads 

E = jcPx Ql A0P + y^gl^P + 1/(101')) (52) 

In the static case jo = e^|(/>p^o so that we can use the Gauss law to 
solve for giving 

^°=2&^ 
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Inserting this in (j52p one can easily accommodate the energy a la 
Bogomol'nyi, a sum of squares plus the a term proportional to the 
magnetic flux, 



E 



j dPx (^{Di±iD2)<p\^ + 



(54) 

where we have discarded a surface term which can be converted to 
a line integral vanishing for finite energy configurations. Then, for a 
fixed value of the flux, <I> = InN/e we have a bound for the energy, 

E>TTe(j)l\N\ (55) 

The bound is saturated when the following first-order equations hold, 

Di(l) = =Fi-D20 

where m = 2e'^(/>Q/|^|. This mass coincides with the vector and the 
scalar field masses m = nij^ = m^, this again signaling some relation 
with N = 2 supersymmetry. 

Being a minimum for the energy, a solution of these equations also 
solve are also static solutions of the eqs. of motion. 

Axially symmetric solutions have been studied using an ansatz 
of the form (I42p . Concerning the Higgs, the field profile is qualita- 
tively analogous tho those for the Maxwell-Higgs system. As for the 
magnetic and electric fields, as it happens for the Maxwell- CS-Higgs 
model, they are concentrated in a ring with its maximum occurring 
when |0p = (/)q. 



Semilocal vortices 

Suppose that one replaces the complex scalar field in the Maxwell- 
Higgs model by an SU{2) doublet cj) = {4>^,4>^^) [25]. The action is 

5 = 1^4^ (-iF^,F^- + ^(Z)^</.^)*(D^0^) + ^(Z)^c/>^^)*(I)V^)- 
V{\4>'\' + \<t>''\')) (57) 
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This action is not only invariant under a local U{l)Local g^iuge group 
but also under a global SU{2)Giobai group. The total symmetry group 
is then 

G = SU{2)giobal X U{l)local (58) 

The potential has an 0(4) symmetry and when it corresponds to a 
symmetry breaking one the vacuum manifold is a three sphere defined 

by 

4' + 4' + 4>i' + 4'' = 4>l (59) 

which is simple connected so that there are no topological string so- 
lutions (ni(S'^) = 0). Now SU{2) global rotations in the direction of 
the constant doublet ((/)q,(/)q^) leave invariant the vacuum so that the 
symmetry breaks to a global U{1)g, 

SU{2)giobal X U{l)local ^ U{l)global (60) 

Flux-tube solutions to the equations of motion of such a model 
have been presented in ^25j where, in view of the symmetry scheme, 
they have been called "semilocal strings" . In order to understand the 
nature of such solutions, let us note that the equations of motion for 
the model are formally the same as those for the Abelian Higgs model 
with just one scalar, except that now such scalar has been replaced 
by an SU(2) doublet. Then any (Nielsen-Olesen) vortex solution 

(j) = f{r) exp(iiVvp) , Ai = a{r)di(t> (61) 

of the ordinary model can be extended to one of the present one by 
introducing an SU{2) doublet $o of unit norm, 

^o=( ^"^'^'K.) (62) 
y V 1 — exp(«99 ) y 

and writing the semilocal solution in the form 

^> = /(r)$o, A = air)d^(P (63) 

Here < a < 1. 

Now, in order to have finite energy, covariant derivatives of and 
should vanish at infinity, this implying 

lim Di4>^ = (idif^ - iNdup) = 

r— >oo 

lim Dicf)^^ = ^ (i9j/^ - iNdiip) = (64) 
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so that the two phases can differ just in one constant 

/ = N^ + C 
ip^^ = Nip 



(65) 



and we have that asymptoticahy, finite energy configurations should 
take the form 



Although there are still a U{l)Giobai rotations of this configuration 
that pass to another one which is still a minimum of the potential, the 
energy of the transformed configuration should be infinity because 
of the change in the covariant derivatives. That is, the space that 
label finite energy configurations does not correspond to the vacuum 
manifold V (the zeroes of the potential) but to the gauge orbit from 
any point $o £ V. In this sense the relevant homotopy group is 
Hi ([/(!)) = Z. Thus configurations with different winding numbers 
are separated by infinite energy barriers, an information that is not 
contained solely in the vacuum manifold V but also in the kinetic 
energy term. Now, since n(V) = it is not topologically guaranteed 
that a configuration characterized by some number is stable. 

In summary, although the vacuum manifold is simply connected 
and a field configuration that winds at infinity can be unwind without 
any cost of potential energy, this does not guarantee that there is no 
cost of the gradient part of the energy and in fact, this is what happens 
due to the non-triviality of Tl{G Local / H Local) = IIi(C/(l)) = Z. 

Stability of the configurations have been studied in [26|-[6]. Again, 
the Bogomol'nyi point plays a central role;for = 8A one can repeat 
the Bogomol'nyi argument and write a bound for the energy, which is 
just that given in eqs. (j25p - (|26p but where the single complex scalar is 
replaced by an SU(2) doublet. Being a local minimum, the solution 
written in (|63p with /(r) and a(r) the radial functions for the Nielsen- 
Olesen solution is classically stable. One should mention the existence 
of zero-modes [56] which will be described below. 

Non Abelian vortices 

The first reference to (local) vortex solutions in non-Abelian gauge 
theories can be found in the Appendix of Nielsen and Olesen pio- 
neer paper |9|. However, it was rapidly recognized that the boundary 




(66) 
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conditions imposed in [9] have solutions that in some cases could be 
deformed continuously to topologically trivial ones [J^. It was then 
understood that in order to have topologically stable solutions, sym- 
metry breaking has to be complete and hence more than one Higgs 
scalar has to be introduced. This fact was discussed in detail in [28j- 
[50] where it was observed that one needs ^i{Gl/H) to be nontrivial. 
For example, for SU{2) gauge theory, just one Higgs field in the ad- 
joint representation is not enough: one needs a second one and a choice 
of symmetry breaking potential so that the only residual symmetry 
group element is just the unit matrix in the adjoint so that H = 7^2 
and ^>i(5C/(2)/Z2) = Z2 (Z^v when the gauge group is SU{N)). 

Let us discuss in some detail the necessary conditions on the vac- 
uum manifold in order to have stable vortices when the gauge group 
is G = SU{N). To begin with, the election of a representation for 
the Higgs field is crucial. Indeed, if we choose the fundamental rep- 
resentation of SU{N) or any other faithful representation (i.e., iso- 
morphic to the group) the fundamental homotopy group Ili{G/H) = 
Ui{SUiN)/I) = = Ui{SU{N)/I) since 5C/(iV) is simply connected. 
The situation is different if the Higgs fields are in the adjoint, which 
does not correspond to a faithful representation of SU{N). 

Noting that elements such that 

n°^ = exj){27rim/N)I , m = 0,1, . . . N - 1 (67) 

can be taken as the elements of the Abelian group Zjy, one can define 
representatives of each homotopy class in the form 

nm{(p + 2tt) = exp{2mm/N)nm{(p) (68) 

In brief, although we say that the symmetry is completely broken, 
there is the residual symmetry H = Z^- and then 

Ui{G/H) = Ui{SU{N)/Zn) = (69) 

It should be noticed that, taken as a gauge transformation group ele- 
ment, Qmi'^) is singular, are the non- Abelian counterpart oiexp{imip), 
for (Abelian) Nielsen-Olesen vortices. 

Let us now discuss the number of Higgs scalars (in the adjoint) that 
one has to introduce to break the symmetry so as to have topologically 
non trivial solutions. We start by recalling the standard Cartan-Weyl 
basis, 

[Hi,Hj] = 0, i = l,2,...,A^-l 
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[Hi, Ea] 


— ^iaEa 


[Ea, Ea] 


— ^iaHi 


[Ea, Ep] 





E^i if r^ = ra + rp 
otherwise 

Here Va is a root vector associated with Ea, ria one of its components. 

It has been proved that in order to completely break the symmetry 
and have topologically non-trivial solution it is enough to have one 
Higgs field ^ in the Cartan subalgebra (the one generated by the 
-ffj's) and — 1 fields in the orthogonal complement |31j-|32j. 
With this in mind, let us take as the gauge Higgs action 

(71) 

where 

= A^t^ , D^ = d^ + e[A^, ] 
E^i/ = d^Ap — dpAp^ + e[A^, Ap\ 

[t'^, t^] = if'bc^c ^ r^^^a^b ^ 1 ^af, ^ a, b = I, 2, . . . , - 1(72) 

The potential should be chosen so as to ensure complete symmetry 
breaking. 

A nonabelian extension of the Nielsen-Olesen ansatz can be then 
proposed in the form 

$^ = r{r)n-\^){Ea^+Ea^)nM , A = 1,2,...,N-1 

N-l 

= - a(r)M , = 

e 

M = diag(l/A^,l/iV,...,l/A^,(l -A^)/A^) (73) 

where M belongs to the Cartan subalgebra. Note that ^ is taken as 
a constant everywhere. One can see that with this ansatz all fields 
satisfy the same equation of motion which, as it happens with 
the gauge field one, coincides with those for the Abelian Higgs model 
so that their properties are the same as those already discussed for 
Nielsen-Olesen vortices. 
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An "electromagnetic" tensor can be introduced in order to 
characterize the vortices. Fohowing [31]- [32] we choose 

_ Tr(MF,.) 

Tr(M2) ^'^^ 

Then, the flux associated to the magnetic field J^i2 reads, for ansatz 

m, 

e = l [ Sxe.^Tij = -^2TTn (75) 
^ J y Ne 



It is important to stress that ansatz (j73p in fact works for any 
n = 2k + m, with m = 1, 2, . . . , — 1 and k £ Z. What is happening 
is that the magnetic flux is not a topological number as in the Abelian 
case. Indeed, as explained above, the topological charge is an element 
of Ztv label in this case by m. 



The Nonabelian Bogomol'nyi bound 

A Bogomol'nyi bound for the nonabelian vortices discussed above was 
found in [33]. We here sketch the SU{2) case which can be easily ex- 
tended to SU{N). As explained above, the potential should be chosen 
so as to completely break the symmetry. One can take a potential of 
the form 

V{^, ^) = ^Tr (^($2 _ ^2^2 ^ gf^^2 _ ^2^2 ^ c,"($^)2) (76) 

minimum for the energy will take place when 

Trim) = (77) 

The Higgs field ^ in (j73p does not play any dynamical role since the 
ansatz corresponds to the condition 

= (78) 

everywhere. Then, the energy E (per unit length) of a general config- 
uration is always bounded by that satisfying the ansatz (j73|) according 
to 

E > <^lTr I (fx (^F^j + A^A^ + +^^(^^ - 1) J (79) 
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Now, as usual, the r.h.s. in (j79p can be written as a sum of perfect 
squares plus additional surface terms/. After some work one finds 

1/- ^^x2 1 ^ ^ 

Di^ + ben"^ A Di^ - — — — -enFii.^ + 



2(1 + a2) '^-^J-J 2(1 + a^^^'' '' 

^ -.2 -,^2 



with a, an arbitrary constant, should be chosen so that the term in 
the last line is semipositive definite. This means 

lal = 1 if A > - 

A if A < ^ (81) 



1 



2(1 + a2) 

Being the last term in the second line of (jSOp related to the magnetic 
flux ([75]) the energy bounded can be written as 

E > ^^vrn if A > - 
- 2 - 8 

E > ^V2X<^l2TTn if A < - (82) 

8 

and the bond saturates whenever the following BPS first order eqs. 
hold 

1 + 

dS + asij^ A Dj"^ = (83) 

Now, one can easily see that these two eqs. are compatible if and only 
if a = ±1 this implying 

A = |! = l (84) 

and we again find the Bogomol'nyi point as a necessary condition for 
bounding the energy. 

In summary, the complete set of BPS equations can be written as 

4 ±^e,,(l-«>2)^ ^ 

Di$±eij^ ADj^ = 

Di^ = (85) 
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with the double sign related to the n sign. At the bound, the energy 
is again proportional to the magnetic flux 

E = $o^|n| (86) 

but, as explained above, \n\ is not the topological charge. Indeed, the 
topological charge is, in this case, an integer modulo 2 and hence the 
more stable solutions other than the trivial vacuum are those with 
n = ±1. 



2 Nonabelian Chern-Simons vortices 

Non-Abelian Chern-Simons vortices have been first discussed in [31]- 
|32j . Let us start by writing the Chern-Simons action for the case of 
an SU{N) gauge theory in d = 3 space-time dimensions, 



Scs[A] = le^'^^Tr j d\ {p^pA^ - '^A^A^A. 



(87) 



This is not a gauge invariant action. Indeed, under a gauge transfor- 
mation 

Al = g-^A^g + '-g'^d^g 
the action changes according to 

Scs[A3]=Scs[A]+87T^IJu;[g] 
with io[g] the winding number associated to the gauge transformation, 

u^ig] = ^Tre"^^ I g-^d^gg-^dpgg-^d^g = me'L (90) 
Since for quantizing the theory one needs that 

eMiScsm) = exp(i5cs[^]) (91) 
(with h= 1) one needs the coefficient jj, to be quantized 

2 

6 771 

^ = , ni G Z (92) 

47r 

Note that this condition is valid both in Minkowski and in Euclidean 
space since the CS term does not depend on the metric and is then 
unaffected by a Wick rotation. 
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We can then add the Chern-Simons term to the Yang- Mills Higgs 
system ()7ip (considered in 2 + 1 dimensions) 

/ N-l 
V ^ A=l 
- V{^,^)j +Scs[A] (93) 

and look for vortex like solutions which, as in the Abelian case, should 
now carry not only magnetic but also electric flux. 

We shall again propose an axially symmetric ansatz as ()73p but 
this time including also the Aq field, 



$^ = r{r)n-\ip){Ea^ + Ea^)nniip) , A = 1,2, . . . , N - 1 



N~l 



= - a(r)M , Ar = 

e 

Tl 

Ao = -ao(r)M 

e 

M = diag(l/Af,l/iV, ...,1/Af, (1 - Af)/A^) (94) 
Given the "electromagnetic tensor" defined as in ()74p 

_ Tr(MF^. ^ 
Tr(M2) 



with a magnetic flux defined as in (|75|) . 



1 /■.2. - 27r 



Q = - d xEijJ^ij = —j=^n (96) 
^ J VNe 

we have now an electric charge Q which can be deflned as follows. We 
first define the electric field from (1951) as 



Ei = Toi (97) 
which can be seen to reduce, with our axially symmetric ansatz to 

Er = T,.. = -^'^ (98) 
ViV dr 
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Then, from the Gauss law, 

dr 

1 



+ a = nB (99) 



where 



B = ^e^JT^J , a = e^fao (100) 

since hmr_>oo -Er = according to the boundary conditions, one gets, 
from (j99p a relation between the charge Q, 

Q = j d^xa (101) 

and the flux, 

Q = (102) 
Now, due to the quantization condition (f92]) . eq. (|102l) becomes 

Q = mnQo , m, n G Z 



Charge quantization can be connected with the angular momentum J 
of the vortex, 

J = J d^xEijXiToj (104) 
with T^jy the energy momentum tensor. One can easily find that 

J = -^ = -^,nm (105) 



Starting from ansatz (|Mj) solutions can be constructed exactly as 
in the previous cases. Both the magnetic field and the Higgs scalar <I> 
have the same qualitative behavior as for neutral vortices. Concerning 
the electric field, it can be related with that of the abelian Chern- 
Simons model already discussed. 

Also as in the Abelian case, one can find Bogomol'nyi equations 
and a bound starting from a Lagrangian where the Yang-Mills term 
is absent and choosing a sixth-order potential for the 4? field |34j . 
restricting to the SU (2) case one starts from 

S = Scs[A] +Tt j d^x {D^^^D^"^" + D^'-^D^,-^ - V{^, ^)) (106) 
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with 

*) = Vi + V2 + c/(*"$")2 (107) 

Again the ^' field will be taken as an spectator just present to com- 
pletely break the symmetry. It will be chosen so as to make vanish V2- 
The third term precisely ensures symmetry breaking. Concerning Vi , 
as in the Abelian case, if one wants to obtain first order Bogomol'nyi 
equations one needs a sixth order potential with a coupling constant 
which is related to the gauge coupling constant 

Vi[<^] = ^^^^^ - $0) (108) 

The energy for a static configuration is given by 

^2 



E = J (fxToo = J Sx {^-{Ao A + \{DSf + ViM 



(109) 



Now the Gauss law, 



e(^ A D°$) = ^EijFij = liB (110) 



allows to write 

^oAI = ^5A$ (111) 
so that the energy becomes 

^ = /^'^ (2^1^!' + V^'^^^ + ^^[i^'O ^^^^^ 

which can be written a la Bogomol'nyi in the form 

{B ^ - T ^$^(^-.5)1 (113) 



We recognize in the last term the magnetic flux 

e = - / (fx4f.B = —n (114) 

so that the energy is again bounded by the flux, 

E > ^lirn (115) 
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and the bound is attained whenever the first order BPS equations 
hold, 



Di^ ± Eij^ A Dj^ = 

^T^I^-l'd^l'-^-o)^ = 

Di^ = (116) 

Quahtatively, the vortex solutions to this system are the same as those 
solving the second-order Euler-Lagrange equations described in the 
precedent section. 

More on non-Abelian vortices and 
semilocal strings 

We have already constructed semilocal vortices for the case in which 
the gauge group was U{l)iocai but one considers an SU {2)giobai. doublet 
of complex scalars. One can extend the construction to a non-Abelian 
gauge group U{N)iocai provided the global group SU{Nf)giobai is such 
that Nf > N. 

Before explaining this, we shall start buy describing a class of non- 
Abelian vortices different from those discussed in sections 5 and 6. 
They were originally discussed in refs.[35]-|37j. Let us consider gauge 
fields taking values in an U{Nc)gauge theory and scalars carrying 
not only a color index a = 1,2, . . . , N but also a flavor index i = 
1,2, .. . ,Nf. We shall first consider = Nf and write these scalars 
as an A^ X A^ matrix, 



^ = m, a,i = 1,2,..., N (117) 



The action is 



S = j (fx (^^TV(F^^F^^) + 'IV(D^$t/)M$) _ ^Tr($tr^$)2 

_^(Tr($t$)_$2)2^ (118) 

Here , ^ = 1, 2, . . . , A^ are the SU (A) generators normalized so 
that Tr(r^r^) = 1 and 



D,^ = {d, - \aI - \aIt^)^ (119) 
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The theory has a U{N)iocai x SU{N)giof,ai symmetry. On the scalars 
this acts according to 



U£U{N)iocal V SUiN) global (120) 

Concerning symmetry breaking, the last term in (jllSp forces <I> to 
develop a vacuum expectation value while the last but one forces it to 
diagonal, 

^vac = ^oSt (121) 

This vev is preserved only for transformations in which U = 
that is, one performs a global gauge transformation and a related 
global flavor transformation simultaneously. The pattern of sponta- 
neous symmetry breaking is then 

U{N)iocal X SU{N)giobal ^ SU{NUag (122) 

This mechanism was known as color-flavor locking in a QCD context, 
[39]-[10]. 

Let us now discuss the vortex solutions of this model. For simplic- 
ity, we shall consider the case N = 2 [36] but the generalization to 
SU{N) will be evident. 

Since there is a U{1) gauge symmetry which is spontaneously bro- 
ken, there are of course Abelian Nielsen-Olesen vortex solutions which 
we have already studied. But there are other novel solutions, related 
to the non- Abelian character of the theory which are obtained as fol- 
lows. 

Let us take the vacuum ()12ip and make just one of its components, 
to wind. This, in the particular N = 2 case means 

^. = ^o(^"f^^ J) (123) 

Once one has such a vacuum, it is easy to write an ansatz for the 
Higgs scalar everywhere, 

^_,„(WOe.pM ^^0^^ ,,,,, 



such that 



lim /i(r) = lim /2(r) = $o 

r— >oo r^oo 

/i(0)=/2(0) = (125) 
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Concerning the ansatz for the gauge field, 



Af(^) = {1 - a{r))d.,ip (126) 

with 

hm a(r) = hm a^ir) = 

r— >oo r—Kxy 

a(0) = aiv(O) = 1 (127) 

This ansatz, in which ^vii is the scalar components that winds is 
usually called a (1,0) string. If one chooses ^v22 to be the sole topo- 
logically nontrivial component will have a (0, 1) string, 

V f2{r)exp{i(p) J 
Let us note that if we define 

ni[<^] = - idi<^) ^-^ (129) 

then, for the (1,0) string one has 

a[^.] = (j l^di^ (130) 

Then, if we associate rij with a topological charge T through the 
formula 

T = Ti^nidx' (131) 

so that 

^(1,0) = ^(0,1) = ^vr (132) 

Let us note that if instead of vacuum (|123|) or (|128|) one proposes an 
"Abelianized" one 

one would obtain for T 

T(i,i) = 47r = T(i,o) + T(o,i) (134) 



25 



Now, since the covariant derivative (jllOp has to vanish at infinity, 



A^L = {d, - - UfT^^ <f.L = (135) 
one easily finds that 

i(a,$.)$-^ = l (^° + AfT^)L (136) 

T = ^Tr ^{Ai + AfT^)dx' = j Aidx' (137) 

Note that only the U{\) component of the gauge field contributes to 
the topological charge We then conclude that the magnetic flux <I> of 
these vortex solutions is precisely connected to T, 

a = -T (138) 
e 

and in accordance, one can see that the Bogomol'nyi bound will turn 
to be 

E > = -^oT (139) 

e 

3 Bogomol'nyi equations and Super- 
symmetry algebra 

We already mentioned at the beginning that the Bogomol'nyi point 

= 8A (140) 

at which solutions to a first-order set of equations (Bogomol'nyi eqs.) 
exist and an exact vortex solution is known, not only corresponds to 
the limiting point separating type-I from type-II superconductivity, 
but, as noticed in [12j, coincides with that which has to be chosen in 
order to have a supersymmetric extension of the model 

This connection between Bogomol'nyi equations and bounds and 
supersymmetry began to be clarified after the work of Witten and 
Olive [H]. Let us then review the results in that paper which starts 
with a two 1 + 1 dimensional supersymmetric model with Lagrangian 

L = -\ ((9^0)' + + V^[<P] + V'[(t>]i;^l^) (141) 
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where tp is a Majorana fermion and V[(J)] an arbitrary function. Our 
metric is diag g^i, = (— 1, 1), (/U = 0, 1). The gamma matrices are 

„o_ „-„2_/0 -1\ i_ , fl 



The charge conjugation matrix C is 

C = a^ , {-ff^fC = -C^i" (143) 

Here a* are the standard Pauh matrices. We write the Majorana 
spinor in the form 

and then 

^ = -^tC = (V-" --0^) (145) 
The equations of motion are, 

^V' - F'V = (146) 
Under the following supersymmetry transformations 
5^ = eip 

5^ = -e{^(j) + V[4>]) (147) 

the Lagrangian changes as a total divergence and then one finds, using 
Noether theorem that the conserved supersymmetry current is 

Ji' = - {^4) + V [({>]) -i^i) (148) 

so that the conserved supercharge is 

Q = - j dx{^(t> + F[<^]) 7V (149) 

Writing 
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one finds from (jl49p 



Q+ = J dx{do4>4^-^ + {di(t> + v[4>])ii^~) 

Q- = j dx (5o(/. V" - - Vm^-^) (151) 
Calling H the Hamiltonian and Z the quantity 



Z = y (152) 

one finds 

{Q^,Q^] = 2{H ±Z) (153) 

where we have used 

{^l}^{x,t),i;^{cy,t)] = 5{x-y) (154) 
Note that Z can be written in the form 

with {dF[(j)])/{d(p) = V[(p]. Then, taking as an example the usual 
symmetry breaking potential with two ground states (p = ita, 



one has 



and so 



V = -{(j)^ - a^) (156) 



Z = X [ dx-^ (a^ct) - ) (158) 



oo 



dx \ 3 



Consider static configurations. This integral vanishes for topologically 
trivial states, 4){x = — oo) = 4){x = oo) = ±a. If instead one takes 
(j){x = — oo) = =F^^ and (j){x = oo) = ita one has Z = ib4Aa^/3. One 
associates the positive value for a kink state and the negative value 
for anti-kink state. Appropriately normalized, the value of T can be 
related to the usual kink topological charge defined as 

T = — / dx— = (j){x = oo) — (j){x = — oo) (159) 
2a dx 
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Eq. (jl53p implies, for a particle with rest mass M = H that 

M > \Z\ (160) 

and the bound is saturated for those states \phys) such that 

Q+\phys) = or Q^\phys) = (161) 

so that physical soliton or anti-soliton states such that (jl62p holds are 
annihilated by (5+ or Q^. Now, eq. ljlGip implies, according to (jl5ip 
(in the static case), 

{di^ + F[(/'])V'- = or {dicj) - V[(l)])Tp+ = (162) 

and then one has the following Bogomol'nyi eqs. 

di(j) = -V[(l)] or 9i(/> = V[(l)] (163) 

Let us note at this point that the energy for static configurations (the 
rest mass M) in the bosonic sector takes the form 

M = j dx^{di(j)f + V^[(I)]) (164) 

and can be trivially written a la Bogomolnyi in the form 

M = J dx^ ((ai0) ± V[^]f T Z (165) 

from which the bound (jl60p can be equally obtained. 

Of course, the solution to these BPS equations also solve the second 
order Euler-Lagrange eqs. Moreover, one can look at eq. (ll63p as the 
square root of the equation of motion for static bosonic configurations. 
One can see this by integrating the equation for configurations such 
that d(j)/dx and V[(j)\ vanish at infinity 

= VV {di(pf = V[(l}f (166) 

Note that since the fermionic supersymmetry transformation law 
()147p can be written, for static configurations, as 

putting this variation to zero is an alternative way to obtain the Bo- 
gomol'nyi eqs without necessity of constructing the supercurrent. 
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One can easily see that there is a zero mode for the Dirac equation 
associated with each Bogomol'nyi equation. These zero modes are 



where Xq constant components. One can understand this by not- 
ing that the kink (antikink) solution of the Bogomol'nyi eqs. (jl63p 
solves, together with ■0 = 0, the coupled Euler-Lagrange eqs. of the 
model. But a supersymmetric transformation of such solution will still 
be a solution and this means that Sip (the SUSY transform of tp = 0) 
will solve the Dirac equation. Such 5^ is precisely (167). 

Let us end by stressing that eqs. (jl53p correspond to an = 2 
extended supersymmetry algebra with central charge Z; the equality 
M = Z when Bogomol'nyi bound is saturated eliminates half of the 
SUSY generators (either Q- or (5-|_) which in turn is associated with 
a Dirac equation zero-mode. 

Supersymmetry and ANO vortices 

We are now ready to see how the vortex BPS equations for the Abelian 
Higgs model can be derived by SUSY considerations. We shall follow 
here refs.[32]-[33]. One has here two possibilities; either one studies 
the N = 2 supersymmetric extension of the d = 2 + 1-dimensional 
model looking at vortices as static solutions of eqs. of motion or one 
analyzes the 2-dimensional (Euclidean) model in which vortices can 
be seen as instanton solutions. Let us start by discussing the d = 2 + 1 
approach. 

The 2 + 1 N = 2 supersymmetric model 

We take as metric diag^^jy = (H )• A two component Majorana 

spinor will be written as 



Spinorial indices are raised and lowered using Ca/s = iSap with e y. = 




(168) 




(169) 



1, 



0" = C^f^^p 



(170) 
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Dirac matrices are chosen as 

7^ = CTi , 7^ = (72 , 7° = (J3 , (171) 

with (7* the standard Pauli matrices. Superspace coordinates are writ- 
ten as with 0" two anticommuting spinor coordinates. 

The Higgs field (f), together with a higgsino tp and an auxihary field 
F are accommodated in a complex scalar superfield, 

9) = (j){x) + ei}{x) - ^eeF{x) (172) 

The gauge field Afj, together with a photino p are accommodated in a 
vector superfield 

v^{x- e) = iA^{x){^>'e)a - eepaix) (173) 

Finally a real scalar superfield will be necessary in order to implement 
the Fayet-Iliopoulos gauge symmetry breaking mechanism, 

S{x; 9) = N{x) + 9xix) - ^99D{x) (174) 

where is a scalar, x is a Major ana fermion and D an auxiliary field. 
Calling Va the (super) covariant derivative, 

'Da = ^+iMad^ (175) 
the spinorial electromagnetic superfield Wa is given by 

Wa = \VpDj:^ = Pa{x)+'-S^uaF'"'{x){-f''e)^-'-9e{^p)a{x) (176) 

With all this the = 1 supersymmctric version of the Abelian Higgs 
model Lagrangian takes, in terms of superfields, the form 

Ln=i = J (fe QvF^T^a - ^ + ^er°) (Pa - ieVa) $ 

—D'^SDaS + V2XS^*^ + r]S^ (177) 

where rj is the real Fayet-Iliopoulos coefficient. Integration out over 9 
one gets the N = 1 supersymmctric Lagrangian in components which 
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takes the form 

Ln=i = -\F^''F^, + ]^d^'Nd^N+]^{D^'(t>rD^cl>+]^D^ 
+\/2AZ?|0p + r}D + + V2XN{F*(j) + Fcj)*) 

76 — — 

+-{:il;p(t>--p^i,(t>*)-V2X{'iiJX(l> + xi^4>*) (178) 

Solving the equations of motion for the auxiliary fields one gets 

L» = \/2A|0p + C, F = -2^/2XN(j) (179) 
Here ^ is the Fayet-Ilopoulos parameter which will be written as 

^ = -^2X4)1 (180) 
so that symmetry breaking takes place for 

\4>\ = 00 > (181) 
With this, Lagrangian ()178p takes the form 

Ln=i = -^F'^"F^, + ia™^iV + ^(Z?'^(/.)*Z?^,/.-4A7V2|^|2 
- 02) +'-p^p+ + W - V2XN{^^|; 

IP _ _ 

+-(V^p0-p^0*)-^/2A(V^X0 + xV'0*) (182) 

This Lagrangian is invariant under the following infinitesimal super- 
symmetry transformations 

'^/o = -^e/..ai^""7''e , 5A^ = -ie^^p (183) 



5N = -ex. 5x = -V2\{\<l,\^ -<l,l)e-ired^N (184) 
<50 = eV^ , 6ip = -z7''eL>^0 - \/8AiV0e (185) 

with e an anticommuting real (i.e. Majorana spinor) parameter. 

It is possible to extend this = 1 supersymmetry to an = 2 
supersymmetry. This can be achieved by considering SUSY transfor- 
mations (|183p - ()185p where the parameter e is replaced by a complex 
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one, ec and real fermions x P ^-re combined into a Dirac fermion 

^ = X-ip (186) 
In terms of this new field the Lagrangian ()182p becomes 

Ln=2 = -^F^"F^, + i5^iVa^iV + ip^(/>rZ)^0-4AiV2|^|2 

-xm^ - 4) + + W - ^^ivv^V- 

-^±^(S,^ + h.c.) + + /..c.) (187) 

where Ti^ is the charge conjugate (the complex conjugate) of S. 

Concerning the SUSY transformations (jl83p - (|185p . they take the 
form, in terms of S, 

SJ: = - Qe^.^F'^^T^ + ^^(101' -<Pl) + i (188) 

5A^ = ^ (e^7^S + h.c.) , 6N = ^{€'E + h.c.) ; , dcj) = e^^P (189) 

5^p = - (iYD^cp + \/8AiV0) e'^ (190) 
Here e'^ is a complex parameter which can be written as 

e'= = eexp(ia) (191) 

and so (jl90p can be interpreted as transformations (|183p - (|185p with 
real e followed by a phase transformation for fermions, 

S exp(ia)S , ip — > exp{ia)ip (192) 

However the last term in p87p is not invariant under this phase rota- 
tion, 

i^T.^cj) + h.c. exp(2ia)'(/;S^(/. + h.c. (193) 

so one needs 

A = IT (194) 



as necessary and sufficient condition for the N = 2 SUSY. This was 
a well-known result holding in general when, starting from an = 1 
supersymmetric gauge model, one attempts to impose a second SUSY 
(see for example ^|): conditions on coupling constants have to be 
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imposed so as to accommodate different = 1 multiplets into an = 
2 multiplet. Eq. (|194p is an example of such conditions. Tliese kind 
of conditions for the Abehan Higgs model were obtained, following 
different routes, in [T3]-[n].[i5]. 

We shall now study the supersymmetry algebra in order to reobtain 
the results above in a way closer to the one used by Olive and Witten 
for their 1+1 model. 

The Noether current associated with supersymmetry invariance 
can be compactly written in the form 

where $ and ^ represent generic bosonic and fermionic fields and is 
included to take into account the possible variation of the Lagrangian 
through a divergence term, 

6L = df.A'' (196) 
From this, the conserved supersymmetric charge Q can be defined as 

Q = [ <fxj'^ (197) 



Writing 

Q = e^Q + Qe^ (198) 

in order to take rid of the infinitesimal Grassman parameter we find 
from Lagrangian (jl87p 

Q = |d2x((-V^F^,7A + ^^iV-|(|</'P-^'))7°S+ 

(iim* - ^A^r) tV) (199) 
Concerning the conjugate charge Q one has 



i;j^^-iJD(t>--N4>)) (200) 

If we restrict the model to the case = and, after computing the 
Poisson brackets we put all fermions to zero we end with the original 
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Abelian Higgs model. Since we are interested in static configurations 
with ^0 = we impose these conditions finding the following anti- 
commutation relation among spinor supercharges, 

{Q«,g^} = 2(7o)i'P° + 5fT (201) 

where 

P' = E = j<fx + \\DM' + jM' - 4?) (202) 
wile the central charge T is given by 

T = j Sx (^^e'^Fi.M'' - <g) + ie'Wi(l){D,(pY) (203) 
with i, j = 1, 2. One can easily see that T can be rewritten in the form 

T = [ (fxd0 (204) 



with 

= e'^ {e^lAj + i(l)*Dj(t>) (205) 

so that using Stokes theorem (and taking into account that Z)j0 ^ 
at infinity) one ends with 

T = e(pl j> Aidx' = e(j)l (^^^ , N G Z (206) 

where the integer n characterizes the homotopy class to which Ai 
belongs. To find the Bogomolnyi bound from the supersymmetry al- 
gebra, let us observe that, being the anticommutators (|2U1|) Hermitian 
one has 

{Qa,Q''}{Q'',Qp}> (207) 
Then, using (j20ip one has 

4£;2 _ > (208) 

or 

E>^-^ (209) 

which is nothing but the bound (|26|) originally obtained by Bogomolnyi 
by completing squares. The bound is attained for states \phys) which 
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are annihilated by the charges. To extract from this condition the 
Bogomolnyi eqs., it will be convenient to write 

Q=(?.^) ' Q = iQ^Q') (210) 



after a little work one finds 

Q+ = -leiJF'^J:+-(^{\cj)\^-v^))J:^-iiD,c^r^P+ + {D2<PriJ- 



2 ' ' V2 

Q- = -^e.jF*^S_-(|(|0|2-z;2))s+-i(Z)i</.)>_ + p20)>+ 

(211) 

Hence, if one combines spinor components so that 

Qi = iQ+ + iQ^) , Qn = {Q+-iQ-) (212) 

one has 

Qi = - {]^e,,r^ + |(|</.|2 - t;2)^ - {i {D^4>r - {D^m ^/ 
Qii = - {]^e,jF^^ - |(|0|2 - ^2)^ s,, _ (, ^D^^Y + (1)20)*) 

(213) 

We see that if the bosonic fields satisfy 

i{Di<i^r -{D2<pr = Q (214) 



Qi vanishes (while Qu not). But if we put ejj = we then have 
Q\phys >= for any configuration satisfying (I214p . This in turn will 
suffice to ensure that the bound (|2U7p is attained. 

Analogously, the super symmetric charge will annihilate physical 
states if 

i (L»i0)* + (1)20)* = (215) 

together with = 0. But ()214p and (j215p are nothing but Bogo- 
mol'nyi equations so that we see that when the Bogomol'nyi bound 
is attained half of the supersymmetry is lost. Again, a fermionic zero 
mode arises associated with this supersymmetry breaking and its form 
can be just inferred by SUSY transforming a trivial (zero) solution as 
in the Witten-Olive example we discussed. 
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Uniqueness of Bogomol'nyi equations 

I will discuss here how the BPS structure of general gauge theories 
(depending on the Maxwell invariant F^'^F^i, and F^^F^jj) coupled 
to Higgs scalars is insensitive to the particular form of the gauge La- 
grangian. Indeed, analyzing their supersymmetric extension, one can 
explicitly understand why this happens. Hence Maxwell, Born-Infeld 
or more complicated non-polynomial Lagrangians all satisfy the same 
Bogomol'nyi equations and bounds which are dictated by the under- 
lying supersymmetry algebra (|46j)- 

Due to the interest aroused by Born-Infeld theories in the context 
of the dynamics of D-branes (see [48j and references there) their Bogo- 
mol'nyi equations for different cases were investigated jSlj (early 
constructions were reported in |52j-[53j). Bogomol'nyi equations were 
found to coincide with those arising in ordinary gauge theories. Con- 
cerning SUSY Born-Infeld theories, they were originally studied in 
[M!-i55j (see also (56]). 

One can quickly understand why BPS equations are not sensitive 
to the dynamics that one chooses for the gauge field by noting that 
they can be derived, in a SUSY framework, by imposing the vanishing 
of (half of) the supersymmetry variations of the gaugino and higgsino 
fields and these variations are formally the same for very different 
Lagrangians. The dynamics associated with the Lagrangian enters 
however through the equation of motion for the auxiliary field D (of 
the gauge field supermultiplet) which appears in the supersymmetric 
transformation law for the gaugino. It is then through D that the 
form of the Lagrangian may in principle determine the form of the 
BPS relations. 

To see this in detail, let us take as an example the Abelian-Higgs 
model in d = 3 dimensions, for which Bogomol'nyi equations were first 
derived [15j.|12j and which we have already discussed in section 2. The 
arguments should hold, however, for other models like for example the 
/SO (3) gauge theory and in other dimensionalities of space-time. 

We shall first consider d = 4 dimensional Minkowski space (with 
signature (-|-, — , — , — )) so that one could also use the results to analyze 
other d = 4 models and then proceed to dimensional reduction to 
d = 3. The gauge vector superfield V is written, in the Wess-Zumino 
gauge, 

V = -Oaf^eAf, + ioeex - mex + ]^e9eeD (216) 

Here is a vector field, A = (Aq) and A = (A") are two-component 
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spinors (a, d = 1, 2) which can be combined to give a four-component 
Majorana fermion and D is an auxihary field. 

From V the chiral superfield can be constructed, 



Wa {y, 9, 9) = -iXa + 9aD-- (a^^a^e)^ F^, + 96 [a'^d^X)^ (217) 

Here A, A, D and F^i, = d^Ay — dyA^ are functions of the variable 
y^i — xf^-\-i6a^9 where is the usual 4- vector position. The SUSY ex- 
tension of (standard) gauge-invariant (Maxwell, Yang-Mills) theories 
are precisely constructed from W by considering W'^ and its hermitian 
conjugate W^. 

Now, as stressed in [SI], another superfield combination enters 
into play if one wishes to construct general gauge invariant SUSY 
Lagrangians. In particular, one needs to consider two superfields X 
and Y defined as 

X = liD^DaW^ + DaD'^W^) (218) 



Y = {D'^DaW^ - DaD'^W^) (219) 

16 

with covariant derivatives given by 

O O O 

Da= ^ + 2i (a''9) —, Da = ^ (220) 

when acting on functions of {y,9,9) and 

C) C) c\ 

Da = —, Da = ^ - 2i (0a^) . ^— (221) 

on functions of {y^,9,9). The only components of these superfields 
having purely bosonic terms are 

W\^ = D^ -^F^'•'F^y-i^F^'^'F^y (222) 

W\^ = D^- ^F^'^F^y + i^Ff^'F^u (223) 

X|o = -{D^-^F''''F^y-iX0X-iX0X) 
X\g§ = i9aP98p{D^ -^F'^'Ff.y -iX0X-iX0X) 
X\em = \9999n{D^ -^Ff"'F^y-iX0X-iX0X) (224) 
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and 

Y\e-e = -^OaPedpi^F''''F^, + X0~X-~X~0X) 

yleeee = leoeeailF^"" F^, + X0X - X0X) (225) 

with F^, = {l/2)e^,^pF'-^. 

A third superfield combination is necessary for constructing general 
gauge invariant SUSY Lagrangian. This combination is W^W"^ with 
its highest component taking the form 

w^w\gg§ = eeee (^{d^ - ^f^uF^^ + {^P^.F^^'^f^ (226) 

Remark that in all three cases (|224p - (j226p . all dependence on the 
curvature F^j^^, and the auxiliary field D is through the combination 

t = ^[D'- If^^'F,.^ (227) 

and this fact will have important consequences in our discussion. Here, 
in order to define a dimensionless variable t we have introduced a 
parameter /3 with the same dimensions as F^i, (i.e. dimensions of a 
mass in d = 4). It corresponds to the absolute field in the Born-Infeld 
theory |57j-|58j as will become clear below. 

As it happens for the last component of W'^W'^, also the term W"^ 
(W^) depends on F^i, and D through the combination (j227p . Indeed, 
the last component in W'^ (W^) contains the term — \F^^'^F^y + 
iF^^Fi"" (D2 _ iF^^F^^ - ^Ff,^F^"') so that the sum of e {0 )integrals 
leads to the well-known SUSY extension of the Maxwell theory. 

We are ready to write a general = 1 supersymmetric Lagrangian 
endowed with gauge- invariance in terms of A, y and W'^W'^ 

irf=4 = J_ I (w^d'e + w^d^e) + \f2arst [d'eiwVyx'Y' 

(228) 

with e the fundamental gauge coupling constant, which has been fac- 
torized in both terms for later convenience. 

The second term accounts for the non-polynomial features of the 
general bosonic theory to be supersymmetrized. As explained in [54], 
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supersymmetry imposes two constraints on coefficients arst- Their ex- 
plicit form will not be relevant for our discussion. What one should 
retain is that expression (j228p gives then the most general Lagrangian 
corresponding to the supersymmetric extension of a general bosonic 
Lagrangian depending on the two algebraic Maxwell invariants F^^F^y 
and Ff'^Ff^y. 

We shall focus now on a d = 3, = 2 supersymmetric theory 
which can be obtained from Lagrangian ()228p by dimensional reduc- 
tion. The standard procedure for dimensional reduction, say in the 
X3 spatial coordinate, implies identifying ^3 with a scalar field N. 
It can be shown that without including a Chern-Simons term, the 
bosonic part of the Lagrangian (|228p can only yield electrically neu- 
tral configurations, so that as long as one looks for self-dual equations 
associated with (static) vortices, the Aq field (as well as the field) 
can be put to zero and so we will do from here on (the case N ^ can 
be equally treated without additional complications). So far, without 
the addition of a Chern-Simons term, no electrically charged vortices 
exist and then the most general gauge field configurations are pure 
magnetic Nielsen- Olesen type soliton solutions. This implies that no 
d = 3 version of the F^^F^'^ functional are available and that we can 
simply identify the field strength with the magnetic field B by 



Once the dimensional reduction is carried on, one ends with the d = 3 
version of the SUSY Lagrangian given in eq. (|228p . As it is well known, 
supersymmetry can be extended from A'^ = 1 to A'^ = 2 in this process. 

From what we have seen, the gauge field dependent terms in the 
bosonic part of this N = 2 supersymmetric Lagrangian can be com- 
pactly written in the form 



-F.yF^"' = 
2 



(229) 



with 



B = ^e,kF^'' i,j 



/• = 1, 2 



(230) 




n=0 



(231) 



where t (defined in (j227p ) now reads 



t = 



D^-B^ 



(232) 
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and c„ are some coefficients which can be computed in terms of the 

^rst 

Concerning the Higgs field sector, in d = 4 dimensions the coupling 
between the scalar Higgs field cp and the gauge field arises from 
the superfield interaction term 

Lf_% = ^eMVW (233) 

where $ is a chiral scalar superfield containing a Higgs field (j), a 
higgsino ip and an auxiliary field F. One can easily see that the part 
of La-cP containing the auxiliary field D is [43] 



Lt-%\D = \d\<P\' (234) 

On the other hand, gauge symmetry breaking can be achieved d la 
Fayet-Iliopoulos so that the complete D dependence of the supersym- 
metric Lagrangian arising from the Higgs coupling to and D is 
given by 

Lt'[A, cb, D] ^ Lf_%\n + iff = ^DM' - e) (235) 

where ^ is a real constant. This Lagrangian remains unchanged after 
dimensional reduction so that we can write the D dependent terms of 
the d = 3 bosonic part of the Lagrangian as 



(236) 



In d = 3 space-time, dimensions of parameters and fields are [(3] = tv?, 
[e] = m2, [^] = m2, [A^] = m, [D] = m? and [0] = ma. Then, for 
dimensional reasons, one can infer that coefficients c'^s can be written 
in the form 

Cn = P^Xn (237) 

where A„ are dimensionless coefficients. 

We can now obtain the equation of motion for D so as to eliminate 
the auxiliary field from the physical spectrum 

E (^-^{D' - B')y-'D +'-m'-e) = (238) 



41 



One can easily see that the only nontrivial solution to eq. (j238p takes 
the form 



B = ±D (239) 

These two equations can be readily combined into one which is nothing 
but the well-honored Bogomol'nyi equation for the magnetic field of 
the Nielsen-Olesen vortices 



B = T^^M'-e) (240) 



This shows that the Bogomol'nyi gauge field equation for vortex con- 
figurations is independent of the particular form of the gauge field 
Lagrangian one chooses since we have proven formula (j240p for the 
general supersymmetric Lagrangian (|228p -l- (l233p . It should be noted 
that the choice of Ai = —1/2, A„ = for n ^ 1 corresponds to the 
usual value of the Maxwell term while the choice Ai = —1/2, A2 = 1/8, 
A3 = 1/32, gives a Dirac-Born-Infeld Lagrangian for the gauge 
field, which, once the auxiliary fields are put on shell takes the form 



Ldbi = -P' (^^1 + ^F.-^F^" - 1 ) (241) 

Let us now analyze the N = 2 supersymmetry transformations 
leaving invariant the three dimensional Born-Infeld SUSY theory. We 
shall not write the complete set of transformations but just those 
which are relevant for the discussion of Bogomol'nyi equations, namely 
those for the higgsino and gaugino (which we call ■0 and S): 

1 ^„.. « , f knF'^ + D A / e , 



5,S = (-e^,„F''V + ^)e 



(243) 

where we call e the Dirac fermion transformation parameter (we have 
already made N = Aq = Q and considered the static case). 

As it is well-known by now, making zero half of the SUSY vari- 
ations associated with the higgsino and gaugino fields, one gets the 
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Bogomol'nyi equations. For instance, by demanding that those gener- 
ated by be zero, one gets the following self-dual equation from the 
higgsino's variation 

S^^i/j = ^ Di(p = iD2(i) (244) 

One should note that this transformation law just depends on the way 
the parallel displacement is defined in terms of the gauge connection 
and not on the explicit form of the gauge field action. One can then 
understand why eq. (1244p is completely independent of the particular 
form the gauge field action is chosen, at least for minimally coupled 
gauge theories 0. Regarding the equation derived from the gaugino 
transformation, 

5,^^ = ]^eijF'^ +D = Q (245) 

it could, in principle, depend on the particular Lagrangian chosen 
through the D term. However, as we have seen (eq. (j239p ). the solution 
to the equation of motion for D takes the same simple form for any 
gauge field Lagrangian since D always enters through the combination 

- B\ 

This feature can be also checked by analyzing the two supercharges 
which can be obtained following the usual Noether construction. As 
it has been shown in [51j for the Born-Infeld case, supercharges Q and 
Q can be always put in the form 

Q = i j cfx^^ n[B,D] {j^B + D) + -^ j d'^x i;^ Pel) 

Q = -ij d^x {B + 7°Z)) n[B, D]^-^-jd^x 7°(^(/.)'f^246) 

with 7i some real functional of D and B which can be computed order 
by order in 1//?^. Furthermore, eqs. (j246p also hold when one considers 
not just SUSY Born-Infeld theory but the general Lagrangian, viz. 
eq. (l228p . Only the actual form of 7i will change, depending on the 
different sets of possible a' s coefficients. What one can easily see is 
that the following formula holds 

oo ^ 

H = Ti-Maxwell + ^ Jm^^'^^IB , D] (247) 

n=l ' 

^FoT an analysis of Bogomol'nyi equations in non- minimally coupled gauge theories, 
see ref.[47]. 
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with 



nn[B,D]\s,^j,, = o 



'^Maxwell — 1 



(248) 
(249) 



It is clear that condition Q\phys) = is satisfied whenever (B + 
7*^-D)e = and ^(pe = 0, independently of the precise form the func- 
tional TC takes. Choosing just the upper component of the transfor- 
mation parameter, e_|_, yields again the two Bogomol'nyi equations 
(|240p . (|244p . Of course, this is to be expected since both {B + 7°L>)e 
and ^(pe, appearing in ()246p . provide the transformation laws of gaug- 
ino and higgsino respectively. 

Concerning the supercharge algebra, when the Bogomol'nyi equa- 
tion B = zizD is used, only the Maxwell part of 7i survives, this show- 
ing again why the BPS structure is not sensitive to the particular form 
of the gauge field Lagrangian. 

In conclusion, we have analyzed the most general Lagrangian cor- 
responding to the supersymmetric extension of a general bosonic La- 
grangian depending on the two algebraic Maxwell invariants F^'^F^u 
and F^^F^j^. This general Lagrangian includes, for a particular choice 
of coefficients, the Born-Infeld supersymmetric Lagrangian, and also 
an infinite class of Lagrangians having causal propagation [5l]. We 
have shown why the Bogomol'nyi relations associated with the bosonic 
sector remain unchanged in spite of the actual form of the gauge 
field Lagrangian: Maxwell, Born-Infeld or more complicated non- 
polynomial Lagrangians all have the same BPS structure. 



Cosmic strings appearing as topological defects in Grand Unified The- 
ory aroused a lot of interest about twenty years ago as a source of pri- 
mordial density perturbations from which galaxies eventually grew. 
In this scenario, since Abrikosov-Nielsen-Olesen strings arise in any 
gauge theory with a broken U{1) symmetry, such strings could appear 
whenever a U{1) symmetry becomes broken as the Universe cools. 

However, after the microwave background results from COBE, 
BOOMERanC and WMAP it was accepted, at the end of the 90's, 
that cosmic strings or other topological defects arising at GUT scales 
could not provide an explanation for the bulk of the density pertur- 
bations. 



4 Including axions 
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The possibility of microscopic fundamental strings in superstring 
theories acting as seeds for galaxy formation was also excluded in 
the context of perturbative string theory [59]: among other problems, 
their tension ^ , close to Planck scale, would produce inhomogeneities 
in the cosmic microwave background far larger than observed. For 
example in perturbative heterotic string theory Gfi = ogut/^^'^ ^ 
10~^ while the isotropy of the cosmic microwave background implied 
Gfi = act/T/ievr < 10-3 [60] 

The question was be recently reconsidered, after the relevance of 
branes and new kind of extended objects was understood (see [61] 
and references therein) opening the possibility to the string tension 
to be much lower, anything between the Planck scale and the weak 
scale. In order to understand the nature and structure of such stringy 
cosmic strings, the embedding of BPS objects in super symmetry and 
supergravity models has become an active area of research so that the 
properties of BPS solitons, their connections with the supersymmetry 
algebra and their cosmological applications have been discussed by 
many authors [33].[62j-[73j. 

Having in mind the study of BPS solitons in a string theory con- 
text, where the axion is almost unavoidable, it is natural to consider 
gauge-Higgs systems in which the axion field is included. In particu- 
lar, an = 1 globally supersymmetric model in d = 3 + 1 dimensions 
consisting of an axion superfield S coupled to WaW°^, with Wa the 
chiral superfield strength, was analyzed in [69] and finite energy cos- 
mic string solutions for the resulting bosonic Maxwell-Higgs action 
coupled to an axion field were constructed. Also, the impact of axions 
on dynamics of a d = 3 + 1 Yang-Mills theory supporting non-Abelian 
strings has been analyzed in [73j. 

We shall consider here this issue taking as an example a d = 2 + 1 
space-time dimensions for which a rich variety of flux tube solutions 
exists already when the axion field is absent. Indeed, as we have 
seen, when gauge fields with dynamics governed by a Chern-Simons 
action are coupled to charged scalars with an appropriate sixth order 
symmetry breaking potential, the model admits BPS equations with 
vortex-like solutions carrying both magnetic flux and electric charge 
[22)-|23j. It should be stressed that in the absence of the Chern Simons 
term, electrically charged vortices with finite energy (per unit length) 
do not exist [74]. Hence, the model we are interested in could show 
novel aspects of charged string like configurations when an axion is 
present, in particular with respect to their application to cosmological 
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problems. Of course, this in the perspective that at high tempera- 
tures, a relativistic four dimensional quantum field theory becomes 
effectively three dimensional. 

The model 

The coupling of an axion to a gauge field with dynamics governed 
by a Chern-Simons action poses some problems |75j-|76j. To discuss 
how can they can be overcome, let us consider the following (2 + 1)- 
dimensional bosonic action, 

S = j d'x [^^f{s)hnD^SF^' + \D,M'' + K"{s)\D^S\'' - W{(p,S)] 

(250) 

where 4> and S = s + ia are complex fields, is a U{1) gauge field, 
and F^^ is defined as, 

F^' = e^'^^d^A^ (251) 

W{(p, (f)* , S, S*) is a potential term, / and K are arbitrary functions of 
s, the real part of S, primes stand for d/ds, k is a constant and 5 is 
a dimensionless parameter (which in the 4 dimensional case is related 
to the Fayet-Iliopoulos term and the Planck mass). Finally, is the 
covariant derivative acting on the fields (j) and S according to 

D^(t) = d^(j) - ieAfj_(j) 

D^S = d^,S + 2i5A^ (252) 

As done in [69j for the 3 + 1 model, we shall identify a in (j250p with 
the axion field and s with a dilaton field. Note that because of the 
definition of the axion covariant derivative, the Chern-Simons term 
appears in action (j250p multiplied by the factor /(s), 

Scs[A,s\ = jj d^xfis)e^'''-A^d,A^ (253) 

The action (j250p is invariant under the local transformation 



s - 


S-2i6A{x) 


s* - 


S* + 2i5A{x) 


- 




r - 




A, - 


A^ + d^A{x) 
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The time component of the gauge field equation of motion is the 
Chern-Simons version of the Gauss law and can be used to solve for 
Ao giving 



An 



kB 



where 



2(e2|,/.|2 + 4(52K") 



1 



B ^= f{s)F,y - e'^{Ai + -d,a)djf{s) 



(255) 



(256) 



The energy can be found from the energy-momentum tensor obtained 
by varying the curved space form of the action with respect to the 
metric, 



(257) 



Integration of the time-time component T^^ gives 
£ = j d^x !^\D,4>\^ + K"\DiS\^ + W{4>,S) + 



4(e2|</>|2 + 4(52 ivT" 



After some work, this expression can be written in the form 
£= ld^x\\D^(l)±iDy^f + K"\D^S±iDyS\^ + 
'e2|0|2 + 4(52i^" 



(258) 



4(e2|0|2 + 452ir" 



B± 



+W- 



± e 



^2/(,)/2 
1 



A6K') 



2| A|2 



+ 4^2^' 



AdK'Y 



±e\ 




A5K') 



+ — SjO 
Id 



djlogfis) 



We thus see that the first two lines in 



(259) 



have been accommodated 



as perfect squares. This, together with an appropriate choice of the 
potential so as to cancel the third line, would lead to a Bogomol'nyi 
bound for the energy given by the magnetic flux ^ appearing in the 
last line, 



d^xF, 



xy 



(260) 
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There is however the term in the forth hne in (j259p with no definite 
sign preventing the obtention of a bound. Only if we put f{s) = 1, 
which corresponds to a normal Chern-Simons action for the gauge 
field (see eq. (j253p ). this term vanishes. In that case one does have a 
bound, 

£ > ±e|(/)oP^' = 27re|0oP|n| (261) 
whenever the potential is chosen as 

W = ^ {e'\^\^ + A6^K") {eM^ - - 46Kf (262) 

The bound is saturated by fields obeying the self-duality equations 

D,S = TiDyS 

K^F,y = ^2 (e^l^p + 45^K") {eM^ - IM^) " ^SK') (263) 

where the upper (lower) sign corresponds to positive (negative) values 
of $ and n. 

As it is well-known, the presence of a Chern-Simons term forces a 
relation between magnetic flux and electric charge [20j. This makes 
the Chern-Simons vortices both magnetically and electrically charged. 
To see this phenomenon in the present case let us write the gauge field 
equation of motion, for the case f{s) = 1 in the form 

^e'^^^F^^ = J>^ (264) 

with J'^ = (/9, J) the conserved matter current and p the electric charge 
density, 

= _2 (e2|<^|2 + 4K"6^) A^' - ie((/>5^0* + (/.*5^(/.) - Addf^a (265) 

We then see that eq. (|255p can be rewritten in the form 

p = -kB (266) 

so that the usual relation between electric charge Q = f d?xp and 
magnetic flux in Chern-Simons theories holds, 

Q = -K^ (267) 

Note that both the Higgs scalar and the axion contribute to the electric 
charge. 
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Supersymmetric extension 

The SUSY extension of the Chern-Simons-Higgs system with a sixth 
order symmetry breaking potential was analyzed in [42] . Let us study 
now the case in which the axion field is also present. 

To do so, we will consider the following d = 3 action, written in 
terms of superfields as 

5susY = -\j (^^(^ + ^^'^^^ (^«^) 

+ 2F($t,$,s,St)) (268) 

Here ;U = 0, 1, 2 and a = 1, 2. Tq, $ and S are spinor, complex scalar 
and axionic superfields and 

= {Da - ieTa)^ 

VaS = DaT. + 2i5Ta 

The lowest component in is the gauge field, that in $ corresponds 
to the Higgs field and that in E is 5 = s + ia with a the axion field. 
Concerning H, K and V , they are functionals of superfields to be 
fixed later. Subindexes in these functionals mean derivatives, thus 
Kgs = dsds*K = d^df^KlQ^Q and so on. They should be chosen so 
that the supersymmetric action (j268p is invariant under the super- 
gauge transformations, 

$ — > e'""^^ 

S — > T.a-26A (269) 

for any real scalar superfield A. Written in the Wess-Zumino gauge, 
i.e., 

Ta\e=o = D'^ra\e=o = 0, 
the spinor superfield Ta is given by 

ra{x,e) = ie'i^nbaAf^ix) - 2e^\a{x) 
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where X{x) is a Majorana spinor, the photino. Then, the spinor field 
strength, defined as 

takes, in terms of component fields, the form 

and satisfies the Bianchi identity, 

D'^Wa = 0. 
The complex scalar super field is defined as 

9) = (l){x) + 9^ipa{x) - 9^F{x) 

where (p stands for the Higgs complex scalar field, ipa is a Dirac 
bispinor, the higgsino and F is a complex auxiliary field. Finally, 
the superfield S which contains the axion S as its lowest component, 
can be reduced to a complex scalar field by exponentiation, 

G = e^ 

The supersymmetric transformations (j269p take F" out from the Wess- 
Zumino gauge. One can however implement a composition of SUSY 
and gauge transformations such that the Wess-Zumino gauge remains 
valid. To do that, the new SUSY transformation for scalar and spino- 
rial super fields are, respectively, 

5^2$ = iT^^Q^^ + iK^ (270) 
J^^S = irf'Qk'L + ik 

where Qa = ida + 9^{'j'^)hadf^, rj"' a Majorana spinor and the real scalar 
superfield K is defined as, 

K = i9''{r)abv'Af, + 9^X''ri,. 

Let us restrict the action ()268p to the case in which 

F(S + S^") = if($1',^>) = 1 (271) 
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Then, written in components, the action takes the simple form 

SsusY = Sb + Sf (272) 
with the pure bosonic action given by 

Sb = J d'x !^^e>'--A^F,^ + \D,,cP\^ + R\D^S\^ - \V^\^ - ^\Vs\^^ 
while the fermionic one is given by 



(273) 



+ 



+ VV' ( -|</'l' -v^^)+xx( —R^ - Vs-, 

26e _ , \ ... /2(5e , \ \R. 



|2 



- ^px[=^<PR + v^,j -rx*i^—'P*R + v^sj-'-^\xx*f 

1 , fRs 4.6^ n\ 1 Viig 

+ rx(^'t>'R-Vf.,']+<h'(^^4>R-VfA\, (274) 

where we have redefined R = Kgs- Here we have taken into account 
that 

1 = j d^xd^e {D"^) Wa 
is a surface term which does not modify the equations of motion 

/ = j d^xd^OD^il^Wa) = j d^xD^D" {J:Wa)\g=o 

and can hence be neglected. 

In order to extend the = 1 supersymmetry to an = 2 super- 
symmetry, we shall follow [42j and allow the transformation parameter 
rj to be complex, i.e., ry* 7^ rj. Now, this is the same as making an 
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A'^ = 1 SUSY transformation followed by a U{1) fermion phase rota- 
tion. Thus, the new transformation for fermions will be 'i/'a — ^ e^"ipa 
and tjj* — >■ e~*"'i/'a ™d the same for Xa and Xa- "^^^ new SUSY trans- 
formations act then as rotations on the fermions and one can then see 
that the only terms which do not respect the extended SUSY invari- 
ance are those on the last three lines in (j274|) . Hence, in order to get 
an A'^ = 2 supersymmetric model we need that, 

where V = V{u, v) = 5 + S*). These equations imply that, 

Vu = -- {eh(u) - 26r{v)) H{u) 

9A 

V^ = {25r{v) - eh{u)) R{v) 

K 

where 

■^r{v) = R{v) , -f-/i(n) = 1. (275) 

av au 

We obviously have h = u — uq = \4>\'^ — \4'o\'^, and since S + S* = 2s, 
then R = K" and r = 2K' where primes stand for derivatives with 
respect to s. From (j273p the potential is, 

which is exactly what we obtained in (|262|) . 

In order to get the Bogomol'nyi bound and the self-dual equations 
one can analyze the supercharge algebra as in Alternatively, 
one can directly consider the component field SUSY transformations 

6a4> = -■0a, ^aS = -Xa, 

Salpb = eabF - i{'y^')abD^(j), 6aXb = eabJ - i{'y^)abDf,S, 

6aF = i{^^tD^i:k + 2\a(t>, SaJ = iij^id^^b + 2Xa, 

(276) 

and their complex conjugated {SX^ = r]*"'5aX^) and reobtain Bogo- 
mol'nyi equations (j263p just by putting all fermion fields to zero in 
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(j276p and then ask the SUSY transformations for ^pa and Xa to vanish 
once the auxiUary fields have been written in terms of dynamical fields 
using their equations of motion. The first condition corresponds to a 
restriction to the bosonic sector, the second one implies that physical 
states are supersymmetry invariant. 



String-like solutions 

We present in this section some vortex solutions to the BPS equations 
of motion. We choose for the Khaler potential the form 

K = -M^log{S + S*) (277) 

As in [22], we shall analyze separately two cases: first, we con- 
sider the case in which the vortex is supported by the Higgs field 
( "^-strings" solutions) , in the sense that at infinity it behaves as in 
the ordinary Nielsen-Olesen vortex, with its winding number linked 
to the magnetic flux. Then, we shall consider the case in which the 
vortex is supported by the axion field, a solution that we shall call an 
"s-string" . In this case it is the axion winding number which is related 
to the magnetic flux. 

In the first case, in order to obtain 0-string solutions we make the 
ansatz [69] 

= /(r)e*"^ S = s(r) - liSmO Ag = (278) 

r 

where n is the topological charge of the Higgs and m is the topological 
charge of the axion. 

It is convenient to work with dimensionless variables by defining 

T = ar, a = 

K 

x{t) = v{t /a) 
y{T) = e6~^s{T/a) 

<r) = 00 V(r/a) (279) 
With this convention the equations read, 

2 

y' = {\m\-\n\x) (281) 

r 

z' = -(l-x)|n| (282) 
r 



53 



where /? = M^/(/)q. From the first two equations we can integrate ?/(t) 
in terms of z{t), obtaining 

y(r) = 2{\n\ - |m|)logT - 21og2;(T) + A; (283) 

with k an arbitrary integration constant. Thus, we end with a system 
of two first-order coupled differential equations for x'{t) and z'{t). 

The boundary conditions for equations (|280p and (|282|) can be 
determined as follows. The function x must vanish at the origin, so 
eq. (1282p implies that z also vanishes at the origin as r'"!. For large 
r, the function x tends to 1, thus z also tends to 1 unless \n\ = \m\. 
In this last case, z — > zq, where zq is the solution of the algebraic 
equation 

4 + 1 ^ = 1 (284) 

° A;-21ogzo ^ ^ 

To solve the differential equations we employ a relaxation method 
for boundary values problem. Such a method determines the solution 
by starting with an initial guess and improving it iteratively. 

Already from ansatz (|278p one can see that the magnetic flux and 
electric charge ar quantized according to 

2tt 27rK 

$ = —n , Q = n 285) 

e e 

There is an interesting property, typical of Chern Simons vortices that 
also holds in our model: both the magnetic field and the (radial) 
electric field are concentrated in rings surrounding the zeroes of the 
Higgs field. 

We also solved the BPS equations when the axion field s tends 
asymptotically to a constants ( "s-string" , [69] ) . In this case the mag- 
netic charge is equal to the topological charge of the axion field, so we 
have as an appropriate ansatz, 

(j) = /(r)e*"^ S = s(r) - 2iSm0 Ag = (286) 

r 

With this ansatz the BPS equations take the form 

|m| V y J \ y J 

y' = --|m|(l-x) (288) 
r 

z' = -{\n\-\m\x) (289) 
r 
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Again, we can integrate y(r) in terms of z{t), obtaining 

y(r) = 2{\n\ - \m\) log r - 2 log z(t) + k (290) 

We see from this equation that for consistency, z(t) ~ r""'" for r — > 
oo in contrast with what happens for the (/)-string. Concerning the 
gauge field boundary condition, one has limT-_>oo a^(T) = 1. 

As a summary of the numerical analysis of the BPS solutions, one 
should note that axionless string solutions found in ([22])- ([23 ) are 
not much modified by the axion which, however, contributes to the 
electric charge of the string configuration. These explicit solutions 
could be of interest in the context of cosmic strings and, due to the 
coupling to the axion and their electric charge, their dynamics could 
be very different of that of ANO vortices. 
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